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ABSTRACT 

We examine radiatively driven mass loss from stars near and above the Eddington limit. Building 
upon the standard CAK theory of driving by scattering in an ensemble of lines with a power-law 
distribution of opacity, we first show that the formal divergence of such line-driven mass loss as a 
star approaches the Eddington limit is actually limited by the "photon tiring" associated with the 
work needed to lift material out of the star's gravitational potential. We also examine such tiring 
in simple continuum-driven models in which a specified outward increase in opacity causes a net 
outward acceleration above the radius where the generalized Eddington parameter exceeds unity. 
When the density at this radius implies a mass loss too close to the tiring limit, the overall result is 
flow stagnation at a finite radius. Since escape of a net steady wind is precluded, such circumstances 
are expected to lead to extensive variability and spatial structure. After briefly reviewing convective 
and other instabilities that also can be expected to lead to extensive structure in the envelope and 
atmosphere of a star near or above the Eddington limit, we investigate how the porosity of such a 
structured medium can reduce the effective coupling between the matter and radiation. Introducing a 
new " porosity-length" formalism, we derive a simple scaling for the reduced effective opacity, and use 
this to derive an associated scaling for the porosity-moderated, continuum-driven mass loss rate from 
stars that formally exceed the Eddington limit. For a simple super-Eddington model with a single 
porosity length that is assumed to be on the order of the gravitational scale height, the overall mass 
loss is similar to that derived in previous porosity models, given roughly by L*/a*c (where L* is the 
stellar luminosity, and c and a* are the speed of light and atmospheric sound speed) . This is much 
higher than is typical of line-driven winds, but is still only a few percent of the tiring limit. To obtain 
still stronger mass loss that approaches observationally inferred values near this limit, we draw upon 
an analogy with the power-law distribution of line-opacity in the standard CAK model of line-driven 
winds, and thereby introduce a power-law-porosity model in which the associated structure has a 
broad range of scales. We show that, for power indices a p < 1, the mass loss rate can be enhanced 
over the single-scale model by a factor that increases with the Eddington parameter as r~ 1+1 / Q i\ For 
lower a p (w 0.5 — 0.6) and/or moderately large T (> 3 — 4), such models lead to mass loss rates that 
approach the photon tiring limit. Together with the ability to drive quite fast outflow speeds (of order 
the surface escape speed), the derived, near-tiring-limited mass loss offer a potential dynamical basis 
to explain the observationally inferred large mass loss and flow speeds of giant outbursts in 77 Carinae 
and other Luminous Blue Variable stars. 

Subject headings: Stars: winds — Stars: early-type — Stars: mass loss 



1. INTRODUCTION 

Massive, hot, luminous stars - those of spectral types 
O, B, and WR - continuously lose mass in strong, radia- 
tively driven stellar winds. For these relatively quiescent 
phases of mass loss (i.e., with mass loss rates ranging up 
to 10~ 5 MQ/yr), the central mechanism for coupling the 
radiation to the outflowing gas is understood to be via 
line opacity, augmented through the systematic Doppler 
shift of the line-scattering by the velocity gradient asso- 
ciated with the flow acceleration. The inherently non- 
linear feedback between the line driving and flow accel- 
eration can be s olved self-consistently via a form alism 
first developed bv lCastor. Abbott, fc Kleinl (^975, here- 



after CAK), with modern extensions providing wind so- 
lutions that are in remarkably good, quantitative agree- 
ment with observational inferences of key wind properties 
like the mass loss rate and wind flow speed, and with the 
scaling of these with stellar parameters like the luminos- 
ity and gravity. 

However, at least some massive stars - observation- 
ally identified as the so-called Luminous Blue Variables 
(LBVs) - appear to undergo one or more phases of much 
stronger mass loss. Perhaps the most extreme example 
is the giant eruption of the massive LBV star 77 Carinae, 
which is estimate d to have cumulatively lost 2-10 be- 
tween 1840-1860 l|Davidson fc Humphrevs1ll997t [Smith 
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12001 . representing a mass loss rate ~ 0.1 — 0.5M Q /yr 
that is about a factor 10, 000 times greater than can be 
readily explained via the line-driven wind formalism. In- 
stead, the observational association of these LBV's as 
being near a n apparent upper limit in luminosity for ob- 
served stars ({Humphreys fc Davidsorll979llT984h has led 
to the general view that the strong mass loss may instead 
stem from the star approaching or exceeding the so-called 
"Eddington limit" , at which even the continuum force 
associated with perhaps just electron scattering exceeds 
the inward force of gravity. 

But a key difficulty in understanding how the approach 
or breach of the Eddington limit might lead to stellar 
mass loss lies in the fact that both the radiative accel- 
eration and gravity have a similar inverse-square scaling 
with radius, implying that their ratio T (the so-called Ed- 
dington parameter) has (in simple 1-D models) a nearly 
spatially constant value throughout the star. As such, 
reaching or exceeding the Eddington limit T > 1 would 
appear to leave the entire star gravitationally unbound, 
and so does not constitute an appropriate description for 
steady-state, surface wind mass loss from an otherwise 
stably bound central star. 

A promising solution to this difficulty lies in relax- 
ing the usual assumption of strictly one-dimensional (1- 
D) spherically symmetric stratification, and considering 
how the lateral structuring - or "porosity" - of a medium 
could lead to an effective reduct ion in the coupling be- 
tween the radiation and matter l)Shavivi ri998. 2003). In 
principal, this can provide a way for quasi-stationary 
wind outflows to be maintained from objects that for- 
mally exceed the Eddington limit. A key insight regards 
the fact that, in a spatially inhomogeneous atmosphere, 
the radiative transport should selectively avoid regions 
of enhanced density in favor of relatively low-density, 
porous channels between them. This stands in contrast 
to the usual picture of simple 1-D, gray-atmosphere mod- 
els, wherein the requirements of radiative equilibrium en- 
sure that the radiative flux must be maintained regard- 
less of the medium's optical thickness. In 2-D or 3-D 
porous media, even a gray opacity can lead to a flux 
avoidance of the most optically thick regions, much as 
in frequency-dependent radiative transfer in a 1-D atmo- 
sphere, where the flux avoids spectral lines or bound-free 
edges that represent spectrally localized regions of non- 
gray enhancement in opacity. In such a case, the force 
opacity can be significantly smaller than the Rosseland 
mean opacity. 

In a gray but porous super-Eddington medium, the 
associated reduction in the coupling of the denser, more 
optically thick regions of the atmosphere can bring the 
effective Eddington parameter below unity, and thus al- 
low for a stably stratified medium. But as the density 
decreases outwards and the individual structures become 
more optically thin, they are again exposed to the full ra- 
diation. In a formally super-Eddington medium, the lo- 
cation where the effective Eddington parameter exceeds 
unity can therefore mark the ini tiation of the n et outward 
acceleration for a stellar wind l)Shavivll2001bj) . 

In principal, the full solution for a wind outflow in 
this case requires self-consistent dynamical solution of 
both the formation of the lateral structure, as well as 
multi-dimensional transport of the radiation to deter- 
mine the effective radiative acceleration throughout the 



thick to thin regions of the medium. As this represents a 
quite challenging and computationally intensive effort, it 
seems appropriate first to seek the development of sim- 
pler, perhaps phenomenological approaches that can pro- 
vide some basic insights into the key properties of wind 
structure and the likely scalings of the resultant mass loss 
rates and flow speeds with relevant parameters. More- 
over, phenomenological approaches have the advantage 
that their results can more easily be implemented in more 
complex systems. 

The aim of this paper is to offer such a simplified, 
heuristic approach, one that takes particular advantage 
of insights and analogies from the much more well- 
established theory for line-driven winds. In the following 
we thus begin (§2) with a general overview of the basic 
properties and scalings of the standard CAK wind for- 
malism for driving by a power-law ensemble of spectral 
lines. We next (§3) discuss how the loss of radiative en- 
ergy, or "photon tiring" , places a fundamental limit on 
both line- and continuum-driven mass loss, and then (§4) 
review how a super-Eddington condition in the deeper 
layers of star can instead lead to convection, pressure in- 
version, and generally a highly structured medium. To 
account for the associated porosity reduction in the ra- 
diative driving, we next (§5) introduce a simple clump- 
absorption formalism that, in conjunction with a new 
"porosity length" ansatz (somewhat analogous to the 
mixing length formalism for convcctivc energy transport) 
allows the estimation of the porosity-moderated mass 
loss rate. Drawing on analogies with the CAK model 
for line-driving, we then (§6) develop a power-law poros- 
ity formalism and apply this to derive scalings for the 
mass loss rate and solutions for the wind velocity law. 
A discussion section (§7) then compares these scalings 
with previous porosity models, and with inferred obser- 
vational properties of the giant eruption in n Carinae, 
including its bipolar form. We conclude (§8) with a gen- 
eral summary and outlook for future work. 

2. RADIATIVELY DRIVEN MASS-LOSS BY LINE AND/OR 
CONTINUUM OPACITY 

2.1. General Equation of Motion 

Consider a steady-state stellar wind outflow in which 
the net acceleration v{dv/dr) in the radial flow speed 
v(r) at radius r results from a radiative acceleration 
g ra d that overcomes the inward gravitational accelera- 
tion, GM»/r 2 , 

dv GM» IdP 

9rad ~ - —, (1) 



dr 



dr 



where the mass density p and the gas pressure P = pa 2 
are related (by the perfect gas law) through the isother- 
mal sound speed a. Using the steady-state equation of 
mass continuity (pvr 2 — constant) to eliminate the den- 
sity p, eqn. Q takes the form 

dv GAf* 



dr 



2a 2 da 2 . , 

g ra d H 1— , (2) 

r dr 



where the terms containing the isothermal sound speed a 
arise from the gas pressure gradient. The square-bracket 
factor on the left-hand-side (LHS) allows for a smooth 
mapping of the wind base onto a hydrostatic atmosphere 
below the sonic point, where v < a. But in radiatively 
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driven winds the pressure terms on the right-hand-side 
(RHS) are generally negligible since, compared to the 
gravitational acceleration term that must be overcome to 
drive a wind, these are of order w s = (a/v esc ) 2 w 0.001, 
where v esc = y/2GM^/R r is the escape speed from the 
stellar surface radius -R*. In the development here, we 
consequently drop these RHS gas pressure terms, but for 
now still retain the LHS factor to allow the option of 
deriving wind solutions that map smoothly onto a hy- 
drostatic atmosphere. 

Since the key to a stellar wind is to overcome gravity, it 
is convenient to define a dimensionless equation of motion 
that scales all the accelerations by gravity, 

(1 -w s /w)w' = -l + T rad , (3) 

where T ra( i = g ra dr 2 /GM*, and the gravitationally 
scaled inertial acceleration is 

r 2 vdv/dr 



w 



GM* 



(4) 



In terms of an inverse radius coordinate x = 1 — i?*/r, 
note that w' = dw/dx, where w = v 2 /v 2 sc represents the 
ratio of wind kinetic energy to the gravitational binding 
v 2 sc /2 = GM»/i?, from the stellar surface radius i?*. 

For isotropic opacity the dimensional radiative acceler- 
ation is set by integration over frequency v of the star's 
radiative flux F v weighted by the opacity k„ 



Qrad 



dv K v F v /c . 



(5) 



In the common, simple case of a continuum dominated 
by free electron scattering, the opacity K e (= 0.2(1 + 
X) cm 2 1 g for fully ionized plasmas with hydrogen mass 
fraction X, which we take here to have the standard cos- 
mic value X « 0.7.) is strictly gray (frequency inde- 
pendent), allowing it to be pulled out of the frequency 
integration. This yields an associated acceleration 



9c 



n e F/c ■ 



47rr 2 c 



(6) 



where L* is the star's bolometric luminosity. In the as- 
sociated gravitationally scaled acceleration, the inverse- 
radius-squared dependence of the flux cancels with that 
from gravity, yielding the so-called Eddington parameter 



AwGM*c ' 



(7) 



which for 1-D radiative transport in the outer envelope 
and surface layers is a spatial constant, set by the ratio 
L*/M* of stellar luminosity to mass. 

This lack of spatial modulation presents a key difficulty 
for driving a steady wind mass loss by electron scatter- 
ing, since a wind requires an outward acceleration that 
transitions from less than gravity in the interior to above 
gravity in the outflowing surface layers. As we discuss be- 
low (§5 & 6), development of lateral structure can lead to 
a porosity moderation of continuum driving that allows 
such a quasi-steady mass loss. 

2.2. CAK Formalism and Scalings for Line-Driven 
Stellar Winds 

The most well-established formalism for radiative driv- 
ing in massive stars invokes the operation of line opacity 



from bound-bound transitions, as first worked out in de- 
tail by CAK. As resonant processes, bound transitions 
have inherently large cross-sections, but only for radia- 
tion within a narrow frequency band of width Ais <C vi 
centered on some line-center resonance frequencies i>x . In 
the dense layers of a star, the diffusion of radiative energy 
means that the flux is reduced in proportion to the en- 
hanced opacity, keeping the effective line force small (see 
eqn. EJ. But within the outflowing wind, the Doppler 
shift of the line-resonance out of the absorption shadow 
of underlying material exposes the line opacity to a less 
attenuated flux, and so can lead to a strong radiative 
force. 

Consider for example the case of a single isolated line of 
resonance frequency v\ within the stellar flux spectrum 
F Vl with integrated strength KiAv/vi = qn e over fre- 
quency width Aja For the simple case of radially stream- 
ing radiation in a supersonic flow, the use of approxima- 
tions introduced by Sobolev (1960) gives the ratio of the 
line acceleration to gravity the scaling 



r„ = W v T e 



1 



te 



(8) 



where W v = viF v j F is a line flux weight, q is the dimen- 
sionless measure of line opacity, and t e = pn e c/ (dv / dr) 
is the radial "Sobolev optical depth" of a line with q = 1 . 
An important aspect of this form is that it allows eval- 
uation of the radiative acceleration in terms of strictly 
local quantities, namely the density /?, velocity gradient 
dv/dr, and line strength q. 

A key advance of the CAK analysis was to introduce 
a formalism for including the cumulative effect of a large 
number of lines of varying strengths q, under the sim- 
plifying assumptions that these are independent. The 
central ansatz is that the flux-weighted number distribu- 
tion of lines is a power law in line-strength 



dN 



1 



5 



dq r(a) 
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(9) 



where T(a) is the complete Gamma function, and a is 
the CAK exponent. Here the line normalization Q is 
related to the usual CAK k parameter through k = 
Q 1 ~ a {vth/c) a /(l — a); it offers the advantages of being 
a dimensionless measure of line-opacity that is indepen- 
dent of the assumed ion thermal speed v t h, with a nearly 
constant characteristic value o f order Q ~ 10 3 for a wide 
range of ionization conditions ijGavlev Ill995j) . 

Integrating the single-line force JHJ over the distribu- 
tion 10 , we obtain for the ratio of the CAK line force to 
gravity 



^CAK = 



C{w>) 



(10) 



(l-a)(Qt e y 

In the latter definition, we have eliminated the density p 
(within t e ) in favor of the mass loss rate M = 4irr 2 pv for 
the assumed steady, spherical expansion, with the line- 
force constant defined by 



C 



1 



1 



Mc 2 



[Qr e 



(ii) 



Note that, for fixed sets of parameters for the star (L», 
A/*, r e ) and line-opacity (a, Q), this constant scales with 
the mass loss rate as C oc l/M". 
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As already noted, the smallness of the dimensionless 
sound-speed parameter w s w 0.001 <C 1 implies that 
gas pressure plays little role in the dynamics of any 
radiatively driven stellar wind. Hence to a good ap- 
proximation, we can obtain accurate solutions by ana- 
lyzing the much simpler limit of vanishing sound speed 
a on ^fwl — -> 0, for which the line-driven-wind equation 
of motion reduces to 



w 



i-r„ 



C{w' 



(12) 



Since the parameters T e and C are spatially constant, 
the solution is independent of radius. For high M and 
so small C there are no solutions, while for small M 
and high C there are two solutions. The CAK critical 
solution corresponds to a maximal mass loss rate, and 
requires a tangential intersection between line-force and 
combined inertia plus gravity, for which 



a C r w„ = 1 



and thus 



with 



< = (i-r e ) 



a 



1 - a 



i-r„ 



l- a 



(13) 
(14) 

(15) 



Using eqn. (|11|) , this then yields the standard CAK scal- 
ing for the mass loss rate 



Mcak = 



1 



_Q£e 

i - r 



(1-q)/q 



(16) 



Moreover, since the scaled equation of motion (|12fl has 
no explicit spatial dependence, the scaled critical acceler- 
ation w' c applies throughout the wind. This can therefore 
be trivially integrated to yield 



w{x) = w(l) x 



(17) 



where w(l) = (1 — T e )a/(1 — a) is the terminal value of 
the scaled flow energy. In terms of dimensional quanti- 
ties, this represents a specific case of the general "beta"- 
velocity-law, 



v(r) 



1 



R* 



(18) 



where here (3 — 1/2, and the wind terminal speed = 
v esc yj a(l — r e )/(l — a) scales with v esc = y/2GM*/ R*, 
the escape speed from the stellar surface radius i?*. 

2.3. Extensions of the Basic CAK Theory 

Modern implementations of CAK theory have in- 
cluded refinements to account for proper integra- 
tion of the finite-angle stellar disk |Friend fc Abbottl 
119861 iPauldrach. Puis, fe Kudritzkil I1986H. a nd ioniza- 
tion s t ate variations with radiu s l|Abbottll980HPauldrachl 
119871 iKudritzki et al. I 11989ft . Together with correc- 
tions for a non-zero ga s pressure (see appendix of 
lOwocki fc ud-Doulal2004|) . these typically only introduce 
order-unity corrections to the above theory, and so the 
scalings derived here still form a good general base to 
compare with continuum-driven formalism we develop 
below. 



Another important extension regards the role of 
multi-line scattering, which occurs whenever the ve- 
locity separation between optically thick l ines o beys 
Av < f no . Analyses bv iFriend fc Castorl l)1982|) and 
iGavlev. Owocki. fc Cranmerl l)1995h show .however, that 
the same CAK scalings still roughly apply as long as the 
spectral distribution of lines is uniform (i.e., follows Pois- 
son statistics) . The total wind momentum then scales as 
MVqq ~ (u oc /At;)i„/c, thus providing a potential line- 
driven model for the dense winds of WR stars, which are 
generally inferred to exceed the single scattering limit 
MUoo = L*/c, sometimes by a factor of 10 or more. 

Finally, because of its relevance to our development 
below of a continuum-driving model, we note in pass- 
ing another modest conceptual modification to classic 
CAK theory, namely to account for the fact that, in any 
discrete distribution of lines, the power-law number ap- 
proximation should be truncated at a line-strength cor- 
responding the strongest, discrete line. Using an expo- 
nential truncation at a line strength q max — Q (at which 
q(dN/dq) w 1), the modified, gravity-scaled CAK line 
force becomes 



lines 



= ^CAK 



(l + l/QU 



^ - (i/Qty- a 



(19) 



In luminous stars with T e within a factor few of unity, 
the scaled CAK line force is of order Tqak ~ 1 — r e , im- 
plying that Qt e « (Qr e /(l-r e )) 1 / Q > 1, since typically 
Q w 10 3 . In such cases, the square bracket correction in 
eqn. (|19|) is therefore very nearly unity. For lower lu- 
minosity stars, however, this square-bracket corrects the 
CAK line-force scaling to impose an overall upper limit 
Times < Qr e . Since a wind requires the driving force to 
exceed gravity, Tu nes > 1, this shows that normal line- 
driven mass loss is only possible for stars with r e > 1/Q. 



"PHOTON TIRING" 



AS A FUNDAMENTAL LIMIT TO 
MASS LOSS 



3.1. Line-Driven Winds near the Eddington Limit 

As a star approaches the classical Eddington limit 
r e — > 1, the standard CAK scalings predict the mass 
loss rate to diverge as M cx 1/(1 - T e ) ( - 1 - a *>/ a , but with 
a vanishing terminal flow speed oc yT — T e . The for- 
mer might appear to provide an explanation for the large 
mass losses inferred in LBV's, but the latter fails to ex- 
plain the moderately high inferred ejection speeds, e.g. 
the 500-800 km/s kinematic exp ansion inferr ed for the 
Homunculus nebula of r\ Carinae (Smith 2003). 

Moreover, of course, such a divergence of the mass loss 
rate is precluded by the finite energy available in the stel- 
lar luminosity L*, which sets a so-called "photon-tiring" 
limit for lifting mass o ut of the gravitational p otential 
from the stellar surface ijOwocki fc Gavlevl ll997). 



M tir = 



La 



L A 



= 0.032- 



Mr. 



R^/M* 



v-J2 GAU/R* ' yr R e /M Q 

(20) 

where Lq = L*/10 6 Lq. Comparison with eqn. i|16fl 
shows that photon tiring would limit CAK winds when- 
ever 



i - r e < Qr e 



-, a/(l-a) 



1 



^esc 



(21) 



For typical parameters Q « 2000 and v^ sc /2c 2 w 10 5 , 
we find that for a = 2/3 photon tiring does not become 
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important until the star is very close to the Eddington 
limit, 

1 - T e < 2000 [2 x 10~ 5 ] 2 = 8 x 10~ 7 ; a = 2/3. (22) 

However, for just a somewhat smaller CAK power index 
a = 1/2, the condition is much less stringent, 

1 - T e < 2000 [10~ 5 ] = 2 x 10~ 2 ; a = 1/2. (23) 

This emphasizes that the CAK mass loss rate is ex- 
tremely sensitive to the power index a, particularly near 
the Eddington limit, where the term within the square 
bracket in eqn. I|16(l has a large numerical value, which 
is then raised to a power that depends on a. 

3.2. Photon Tiring in Continuum-Driven Mass Loss 

Let us now examine photon tiring for continuum- 
driven flows, ignoring the effect of line-opacity, and sim- 
ply considering that the continuum opacity now has 
a known radial dependence n c (r). Specifically, we as- 
sume this opacity increases outward until some layer 
becomes super-Eddington, i.e., T c > 1, where now 
r c (r) = K c (r)L*/47rGM*c is the generalized continuum 
Eddington parameter. Defining for convenience the ra- 
dius r — i?» to be the radius at which r c (i?») = 1, 
let us assume moreover that this represents the sonic 
point of an initiated steady-wind outflow. The density 
p» and sound speed a* at this point set the mass loss rate 
M = 4irRlp*a*, but otherwise gas pressure terms have 
negligible effect in the further supersonic acceleration of 
the outflow. 

In terms of the scaled wind energy w and scaled inverse 
radius x defined in §2.1, eqn. again represents the 
scaled equation of motion, with T rac i = r c (x). Note, 
however, that in this form in which T c (x) is presumed to 
be a known spatial function, the mass loss rate itself has 
scaled out, so that the resulting velocity law would be 
entirely independent of the amount of mass accelerated. 

More realistically, as noted above, a given radiative 
luminosity can only accelerate a limited mass loss rate 
before the energy expended in accelerating the outflow 
against gravity would necessarily come at the expense of 
a notable reduction in the radiative energy flux itself. To 
take account of this photon tiring, we must reduce the 
radiative luminosity according to the gained kinetic and 
potential energy of the flow, 

<P_ GM* _ GM* 
2 + R* r 



L{r) = L» 



M 



(24) 



which implies for the tiring-corrected Eddington param- 
eter that 



Frad(x) = r c (x) [1 - m(w + X)} , 

where the gravitational "tiring number" , 



M 



MGM* 



0.012 



M_4V^ 000 



(25) 



(26) 



characterizes the fraction of radiative energy lost in 
lifting the wind out of the stellar gravitational poten- 
tial from .R*. The last expression allows easy eval- 
uation of the likely importance of photon tiring for 
characteristic scalings, where M_4 = Af/10~ 4 M Q /yr, 
Lq = L»/10 6 L©, and Viooo = v esc /1000 km/ s rj 
0.62(M*/R*)/(M Q /R Q ). 
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Fig. 1. — a. Wind energy w vs. scaled inverse radius x(= 
1 — R-tJr), plotted for Eddington parameter r c (x) = 1 + yfx with 
various photon tiring numbers m. b. Same as (a), except for weak 
tiring limit m <^ I, and for various constants c in the Eddington 
parameter scaling T c (x) = 1 + \fx — 2cx. 



Applying the tiring-corrected Eddington parameter 
eqn. I|25|l into the dimensionless equation of motion |J3}, 
we find 

1 - — ) ^ = -1 + T c (x)[l - m{w + x)}, (27) 
w / dx 

where for typical hot-star atmospheres the sonic-point 
boundary value is very small, w(0) = w s = w* = 
a%R*/2GM* < 10~ 3 . Through most of the wind ac- 
celeration, i.e., except near the sonic point itself, the 
sound-speed term with w* can thus be neglected. In 
the idealized limit w* — > 0, we can use integrating fac- 
tors to obtain an explicit solution to w(x) in terms of the 
integral quantity T c (x) = J Q dx'T c (x'), 

1 



w(x) 



m 



1 



-mT a (x) 



(28) 



As a simple example, consider the case 1 with T c (x) = 
1+y/x, for which f c = x+2x 3 ^ 2 /3. Fig. la plots solutions 
w[x) vs. x from eqn. I|28|) with various m. For low m, the 
flow reaches a finite speed at large radii (x = 1), but for 
high m, it curves back, stopping at some finite stagnation 

1 The choice of these functions is arbitrary, to illustrate the 
photon-tiring effect within a simple model. More physically moti- 
vated models based on a medium's porosity are presented in § 6.4. 
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point x s , where w(x s ) = 0. The latter solutions represent 
flows for which the mass loss rate is too high for the given 
stellar luminosity to be able to lift the material to full 
escape at large radii. By considering the critical case 
w(x = 1) = 0, we can define a maximum mass loss rate 
from m max , given from eqn. (9) by the transcendental 
relation, 



= 1 



»r (l) « i _ e 2-2f c (i) 



(29) 



where the last expression provides a simple explicit ap- 
proximation 2 for any realistic r c (l) > 1. Note that re- 
gardless of how large T c (l) becomes, it is always true 
that m max < 1, simply reflecting the fact that the mass 
loss is always limited by the rate at which the radiative 
luminosity can lift material out of the gravitational po- 
tential from R*. For a given photon tiring number m < 1, 
the minimum integrated Eddington parameter required 
to ensure full escape from the potential is 

= . , — ln(l — m) 

r c i = — ^ '- . (30) 

m 

Even without photon tiring, a similar stagnation can 
occur from an outward reduction in the radiative driving. 
In the limit of negligible tiring m <C 1, the flow solution 
<|28l) simplifies to 

w[x) i=s r c (a;) — x. (31) 

For a limited super-Eddington domain, the critical case 
of marginal escape with zero terminal velocity, w(l) = 0, 
is now set in general by the condition f c (l) = 1. For 
example, consider the specific case of a nonmonotonic 
r c (x) = 1 + \fx — 2c x, for which then T c (x) = x + 



2_x 3 / 2 /3- 



Fig. lb plots results for various c. For all 



r c (l) < 1 (i.e., c > 2/3), the material stagnates at the 
radius where T c (x s ) = x s = 4 /9c 2 , and so cannot escape 
the system in a steady-state flow. In a time-dependent 
model, such material can be expected to accumulate at 
this stagnation radius, and possibly eventually fall back 
to the star. This represents another way in which, in- 
stead of a steady outflow, a limited super-Eddington re- 
gion could give rise to an extended envelope with either 
a mass circulation or a density inversion. 

4. STELLAR ENVELOPE CONSEQUENCES OF 
BREACHING THE EDDINGTON LIMIT 

4.1. Convective Instability of Deep Interior 

It should be emphasized that locally exceeding the 
Eddington limit need not necessarily lead to ini- 
tiation of a mass outflow. As first shown by 
I.Toss. Salneter. fc Ostrikerl l|1973j) . in the stellar envelope 
allowing the Eddington parameter T — > 1 generally im- 
plies through the Schwarzschild criterion that material 
becomes convectively unstable. Since convection in such 
deep layers is highly efficient, the radiative luminosity is 
reduced, thereby lowering the associated radiative Ed- 
dington factor away from unity. 

This suggests that a radiatively driven outflow should 
only be initiated outside the region where convection is 
efficient. An upper bound to the convective energy flux 
is set by 



. I dU/dr <aH dP/dr 



a 3 p, 



(32) 



2 This has a maximum error of ca. 
plex, explicit approximation is m max 



which is accurate to within a half percent. 



i. A somewhat more com- 
1 - exp[5(l - r c (l) 2 / 5 )], 



where v conv , I, and U are the convective velocity, mixing 
length, and internal energy density, and a, H, P, and 
p are the sound speed, pressure scale height, pressure, 
and mass density. Setting this maximum convective flux 
equal to the total stellar energy flux L^/Airr 2 yields an 
estimate for the maximum mass loss rate that can be 
initiated by radiative driving, 



M < = .} I, 
a 



max,conv 



2a? 



'-M t , 



(33) 



where the last equality emphasizes that, for the usual 
case of a sound speed much smaller than the local escape 
speed, a <C v esc , such a mass loss would generally be 
well in excess of the photon-tiring limit set by the energy 
available to lift the material out of the star's gravitational 
potential (see ean. I2UI) . In other words, if a wind were 
to originate from where convection becomes inefficient, 
the mass loss would be so large that it would use all the 
available luminosity to accelerate out of the gravitational 
potential. 

A central conclusion here is accordingly that, while 
convective transport provides an alternative to a super- 
Eddington condition in the deep interior, it cannot be 
the regulation mechanism that would allow for a smooth 
transition to a steady wind mass loss of the near-surface 
layers. 

4.2. Hydrostatic Pressure Inversion in a 
Super-Eddington Layer 

Even above the inefficient-convection radius, a limited 
super-Eddington domain could, instead of an o utflow, 
merely induce a pressure inversion layer (Maeder 1989), 
set by integrating the equation of hydrostatic equilil- 
brium (cf. eqn. ^) 

d\nP r-i 

where H = a 2 r 2 /GM* is the usual gravitational scale 
height. As an example, for a narrow (Ar <C r), isother- 
mal, super-Eddington layer, the pressure would increase 
by a factor exp[(Ar/_ff)(r — 1)], where T is the average of 
r over the layer. This exponential pressure increase im- 
plies, however, that such inversions are only possible over 
a limited domain, since eventually the star must match 
an outer boundary condition of negligible pressure. 

Since more realistically the temperature should be ex- 
pected to decline outward, such a pressure inversion 
would imply an even stronger outward increase in den- 
sity. If the super-Eddington condition persists, then per- 
haps another switch to convective transport could again 
reduce the radiative flux, and bring the Eddington pa- 
rameter back below unity. But again this should become 
inefficient at a layer that cannot maintain an outflow 
against photon tiring. This implies that any such outflow 
initiated from the region that convection becomes inef- 
ficient would necessarily stagnate at some finite radius. 
One can imagine that the subsequent infall of material 
would likely form a complex spatial pattern, consisting 
of a mixture of both downdrafts and upflows, perhaps 
even resembling the 3-D cells of thermally driven con- 
vection. Overall, we thus see that a star that exceeds 
the Eddington limit is likely to develop a complex spa- 
tial structure, whether due to local instability to con- 
vection, or to global instability of flow stagnation. Ir- 
respectively, an inversion or efficient convection cannot 
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offer a solution giving a photospheric luminosity which 
is super-Eddington. 

4.3. Lateral Instability of Thomson Atmosphere 

Dating back to work bv lSpiegell Ijl976lll977|) there have 
been speculations that an atmosphere supported by ra- 
diation pressure would likely exhibit instabilities not un- 
like Rayleigh- Taylor, associated with the support of a 
heavy fluid by a lighter one, leading to formation of 
"photo n bubbles" . R e cent quantita tive stability a nal- 
yses bv lSpiegel fc~Taol (fl999^ and bv lShavivl l|2001a|) do 
lead to the conclusion that even a simple case of a pure 
"Thomson atmosphere" - i.e., supported by Thomson 
scattering of radiation by free electron - would be sub- 
ject to intrinsic instabilities for development of lateral in- 
homogeneities. The analysis bv iShavivl l)2001al) suggests 
in particular that these instabilities share many similar 
prop erties to the excitation of strange mod e pulsations 
(e.g.. lGlatze!"lll994l: IPanaloizou et allll997fl . For exam- 
ple, they are favored when radiation pressure dominates 
over gas pressure. Both arise when the temperature per- 
turbation term in the effective equation of state for the 
gas becomes non-local. In strange mode instabilities, the 
term arises because the temperature in the diffusion limit 
depends on the radial gradient of the opacity perturba- 
tions. In the lateral instability, the term depends on the 
lateral radiative flux which arises from non-radial struc- 
ture on a scale of the vertical scale height. 

Note that when conditions of a pure Thomson atmo- 
sphere are alleviated, even more instabilities arise. There 
are of course the aforementioned strange mode instabil- 
ities, which require a non-Thomson opacity. If mag- 
netic fields are introduced, even more instabilities can 
plav a role <IArons)ll992l iGammie 1119981 iBegelmanlfeOOa 
iBlaes &: Socrates! 20031) . We stress however that the 
physical origin of the instabilities is not important to 
our discussion here. The essential point is merely that 
as atmospheres approach the Eddington limit, non-radial 
instabilities do exist to make the atmospheres inhomoge- 
neous, while the typical length scale expected is that of 
the vertical scale height. 

5. SUPER-EDDINGTON OUTFLOW MODERATED BY 
POROUS OPACITY 



Shaviv ( 2000. 2001bl) has applied these notions of a lat- 
erally inhomogeneous radiatively supported atmosphere 
to suggest a new paradigm for how quasi-stationary 
wind outflows could be maintained from objects that 
formally exceed the Eddington limit. A key point re- 
gards the fact that, in a spatially inhomogeneous at- 
mosphere, the radiative transport will selectively avoid 
regions of enhanced density in favor of relatively low- 
density, "porous" channels between them. This stands 
in contrast to the usual picture of simple 1-D, gray- 
atmosphere models, wherein the requirements radiative 
equilibrium ensure that the radiative flux must be main- 
tained independent of the medium's optical thickness. In 
2-D or 3-D porous media, even a gray opacity should lead 
to a flux avoidance of the most optically thick regions, 
much as in frequency-dependent radiative transfer in 1- 
D atmosphere, wherein the flux avoids spectral lines or 
bound-free edges that represent localized spectral regions 
of non-gray enhancement in opacity. The associated re- 
duction in the effective opacity might therefore provide 



a mechanism for the transition from an effectively sub- 
Eddington to super-Eddington condition, and perhaps 
thereby allow an appropriate regulation for steady wind 
mass loss of the outer layers. 

5.1. The Porosity Length and Effective Opacity of a 
Clumped Medium 

The radiative transport in such a complex, 3-D 
medium is likely to be extremely complicated, but to 
model the basic elements of this porosity effect, we can 
consider a simplified picture of a medium consisting of 
an ensemble of localized clumps or blobs. 3 Assuming 
for now that the blobs all have the same characteristic 
length I and mass mj, the characteristic blob density is 
Pb ~ mb/l 3 , implying in a medium with opacity k a char- 
acteristic blob optical depth w Kpbl- If the blobs are 
optically thick, rj, ^> 1, then their effective cross section 
to impingent radiation is just a c g « I 2 . But more gen- 
erally, since for arbitrary optical thickness the fraction 
of impingent radiation attenuated by each blob should 
scale as 1 — exp(— r&), this effective cross section can be 
written as a c g ~ l 2 [l — exp(— t&)]. From this, we can 
accordingly define an effective opacity of the blobs as 

K e g = ~ — [l-exp(-Tfe)] w K- — exp ( Tb , (35) 
m,b nib n 

In the limit that the blobs are optically thin, T& <C 1, this 
effective opacity recovers the microscopic value, K c ff ~ 
k(1 — Tb/2) — > k. However for optically thick blobs, Tb ^S> 
1, the opacity is effectively reduced by a factor l/r&, i.e. 
«eff — * n/n ~ l 2 /m b . 

Let us thus consider a medium that consists entirely of 
an ensemble of such blobs, with a characteristic separa- 
tion scale L I. Then the mean density of the medium 
is given by p « mb/L 3 = (l/L) 3 pb, and the blob optical 
thickness can be written as the ratio 



P 

n = — 

Pc 



(36) 



where the critical medium density at which the blobs 
have unit optical depth is given by 



(37) 



I 2 1 
^ c kL 3 nh ' 
with the latter equality defining a characteristic "porosity 
length" h = L 3 /l 2 . 

As developed further below, this porosity length turns 
out to be a key parameter for determining the nature and 
consequences of the porosity in a structured medium. 
The above shows already that it can have a variety of 
interpretations: 

• As defined here, it is the ratio of the volume per 
blob to the projected surface area of the blob, h — 
L 3 /l 2 . 

• Equivalently, it is the blob size divided by its vol- 
ume filling factor, h = 1/(1/ L) 3 . 



3 For comparison. IShavivl 1200 lbl) considered two limiting cases, 
of either blobs in the optically thin limit or thick but vertically 
elongated, and a third case of a two phase Markovian mixture. Our 
analysis here generalizes the first two cases, but is rather different 
than the statistical description in the third case. Nonetheless, the 
comparison with the third case reaffirms that the mass loss scaling 
does not depend on the exact assumed geometry. 
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• As noted above, the inverse of the porosity length 
times the opacity defines a critical mean density 
at which the blobs become optically thick, p c = 
X/{hK). 

• For mean densities above this critical density p 3> 
p c , it represents the photon mean- free- path of the 
porous medium, h = l/(n c sp) = l/( K Pc), which 
consequently becomes independent of this mean 
density. 

5.2. Mass Loss Rate for Model with a Single Porosity 

Length 

Let us now apply this simple picture of a porous 
medium to model the porosity-moderated mass loss of 
a super-Eddington atmosphere. The generic mass loss 
rate can be written as 



M = 4nR;p*a 
= 4:TTRlp*a 



1 



HGM* 



p c nh 
h J If J a*c 



Rial 



TA-kGM*c 



(38) 



where the basic definition in the first line is multiplied in 
the second line by a series of unity factors (in square 
brackets), which are defined in terms of the porosity 
length h, the critical density p c (= 1/nh), the gravita- 
tional scale height H (= a*Rl/GM*), and the Eddington 
parameter T (= kL*/ ' AttGM^c). After several cancella- 
tions, the result in the third line yields three dimension- 
less scale factors (in parenthesis) that multiply the di- 
mensional scaling set by the luminosity divided by the 
product of the speed of light times sound speed. For fu- 
ture reference, we refer to this last factor as the "basal" 
mass loss rate, with notation 



a*c 



(39) 



So far this is just a simple recasting of the generic mass 
loss scaling, but it now allows us to apply properties of 
a specific porosity model to derive associated physical 
scalings of mass loss from a porosity-moderated, super- 
Eddington medium. 

Let us then first define a scaling for the porosity length 
h. For this, consider arguments analogous to those tra- 
ditionally given for representing convective energy trans- 
port in terms of a characteristic "mixing length" , which 
is generally assumed to scale in proportion to the grav- 
itational scale height H . Specifically, let us make an 
analogous ansatz that the porosity length should like- 
wise scale with this pressure scale height. We thereby 
define a dimensionless 11 porosity-length parameter 1 ' , 



(40) 



which is analogous to the mixing-length parameter, a = 
Imix/ H , traditionally assumed for treating convective en- 
ergy transport. 

Note that this scaling relates the porosity length to 
the pure gas pressure scale height, ignoring any modifi- 
cation due to radiation pressure. A basic justification 
for this lies in the expectation that the effective Ed- 
dington parameter will be substantially reduced in the 



porous medium, so that in the final self-consistent state 
the medium again becomes stratified on a scale propor- 
tional to H. More generally, any effect of radiation pres- 
sure in setting this equilibrium can be accounted for by 
taking this porosity-length parameter to be a function of 
the Eddington parameter, r][T]. However, for simplicity, 
we typically take 77 to be a fixed independent parameter 
in the analysis here. 

Let us next derive the scaling for the sonic density 
ratio, p*/ p c - Given a value of the Eddington parameter 
r > 1 defined from a microscopic continuum opacity k, 
the reduced, effective Eddington parameter for material 
with arbitrary density p is simply given from application 
of eqns. and i|36fl. 



r eff = ^£r« ^rd-e-'/" 

K p 



(41) 



The reduction at large densities now allows a base hy- 
drostatic region where r e ff < 1. The transonic wind is 
then initiated at the point where r c ff = 1, which when 
applied to eqn. (|41|l defines an implicit relation for the 
sonic point density . An explicit solution is given by 



Pc 



r + Wo [-r exp(-r)] « r - l/r (42) 



where Wo represents the princip al branch of the Product- 
Log ( a.k.a. Lambert) function 4 (Jef frey. Hare, fc Corliss! 
1996), and the latter, much simpler, approximate form 
is valid within 6% of the exact solution for all the T > 1 
that are of interest. 
Application of the approximate form from l|42|) into 
yields the mass loss scaling, 



M, 



por ' 



1 \ L 



1 \ L» — L Edd 



W[T] 



— LEdd 

a*c 



(43) 
(44) 
(45) 



The second equality recasts this scaling in terms of a 
difference of the luminosity from the Eddington lumi- 
nosity, L Edd = Att GM^c/ k. This is quite similar to the 
scaling derived bv iShavivl l)2001o|) . with the final equal- 
ity introducing his "wind function" (cf. his eqn. 30). In 
the present formalism this is seen to have the approx- 
imate scaling W(r) w (r + l)/r?7, which hence varies 
from W ~ 2/rj for the marginally super-Eddington case 
r — 1 <C 1 to W ~ l/r) for the strong super-Eddington 
limit T > 1. 

In terms of the luminosity-proportional form (|43|l , the 
dimensional values of the mass loss rate for this singlc- 
porosity-lcngth model scale with the basal rate defined 
in eqn. i!39l 



r]M p 



i - i/r 2 



L* 
a*c 

= 0.001 



M G) L 6 

— — , 46 

yr a 20 

where 020 = a*/20 km/s and L$ = L*/1O 6 L0. Dividing 
this basal mass loss rate by the tiring-limited value in 

4 http : / /mathworld . wolfram . com/LambertW-Function . html 
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(47) 



eqn. i|20[l gives for the associated "tiring number" 
M* v 2 esc 0.032 M,/R, 

ffl = — = i±i- = . 

Mur 2a * c fl 20 Mq/Rq 

We thus see that, under the canonical assumption that 
i) w 1 (i.e., h w H), this single-porosity- length model 
yields a maximum mass loss rate that is a few percent of 
the photon tiring limit. 

5.3. Two- Component Model of Clumps in a Smooth 
Background Medium 

As a next step in developing this phenomenologi- 
cal model for porosity, consider the more realistic case 
in which these localized blobs are embedded within a 
smooth, unclumped component. Let us again assume the 
clumped component consists of many individual blobs 
with mass m b confined to a small size I <C L com- 
pared to the interblob spacing L, giving them a density 
Pb = mb/l 3 that is much larger than the interblob den- 
sity pi — rrii/L . The density averaged over blob and 
interblob medium is 



P 



in 



= Pi+p b (l/Lf, 



(48) 



where the total mass m = mi + m b . 

Again taking the medium to have a fixed microscopic 
opacity k, the individual blobs thus now have an optical 
thickness 

Tb = npbl = Kmb/l 2 = KpfL 3 /l 2 ~ Kpfh, (49) 

where / = rrib/rn is the blob mass fraction, and h = 
L 3 /I 2 is again the characteristic porosity length. As long 
as the blobs are optically thin, < 1, both the blobs 
and the interblob medium have this same effective mass- 
absorption coefficient, namely k. But if the blobs become 
optically thick, T& > 1, then, compared to a uniform 
medium, the effective opacity of this clumped medium is 
now reduced by a factor 



k = 



Kcff 
K 

Km; 



K b m b 



1-/ + / 



l 



n 



(50) 

Here k m i n = 1 — / represesent a minimum value of 
this effective opacity ratio associated with the interblob 
medium, while k b represents the reduced opacity of the 
blob component. 

It is again instructive to consider a few limiting cases. 
First, for optically thin blobs n -C I , we do indeed re- 
cover that there is no effective change 



I - 



fn 

2 



+ 0(t 2 )^1 ; r b «I. 



(51) 



In the opposite extreme that the medium consists almost 
entirely of optically thick blobs, with very little interblob 
material, 1 — / <C I/tt,, we again find that effective ab- 
sorption is reduced by a factor given by the inverse of 
the individual blob optical thickness, 

k -> - ; 1/(1 - /) » n » f, (52) 
n 



which in principal can imply quite a substantial reduc- 
tion in absorption. In the somewhat intermediate case of 
optically thick blobs with an interblob medium that has 
a non- negligible mass fraction 1 — / ^> l/rfc, the absorp- 
tion is dominated by this interblob material, yielding 



fc 



n » 1/(1 - /) > 1. 



(53) 



This last case represents an important component of 
added realism for the model, since it implies that, un- 
like the pure-clumped case considered above, the poros- 
ity reduction does not continue to scale in proportion to 
the inverse density to arbitrarily small values, but rather 
saturates to a minimum or "floor", k m i n = 1 — /. More 
realistically, instead of identical clumps, it seems likely 
that the structured component may be better modeled 
as consisting of a broad distribution in clump properties, 
as is developed in the next section. 

The general picture then is that, in such a porous 
medium, the radiative flux, and thus the radiative ac- 
celeration, should be weighted towards the lower-density 
regions, which, in principle can thus have an even greater 
(super-Eddington) outward driving. However, a further 
fundamental assumption in the analysis in this paper is 
that overall medium can still be described in terms of a 
single-fluid model. One justification lies in the relatively 
high overall densities in the regions of flow initiation - 
which thus implies a strong gas collisional coupling to 
share momentum between lower and higher density com- 
ponents. Moreover, since many of the models for struc- 
ture formation involve travelling wave solutions, the ac- 
tual material can even alternate in time between dense 
and rarefied regions, thus implying a further overall aver- 
aging in radiative momentum addition. With this over- 
all assumption of strong effective momentum coupling to 
maintain a single-fluid medium, we thus leave to future 
studies the issue of possible dynamical differentiation be- 
tween regions of lower and higher density. 



6. SUPER-EDDINGTON WIND MODERATED BY A 
POWER-LAW POROSITY 

6.1. A Power-Law-Porosity Ansatz 

In considering options for further development of this 
basic formalism, let us first note here an interesting sim- 
ilarity in how both continuum porosity and line opacity 
lead to a reduced effectiveness in absorption as a function 
of the relevant optical thickness parameter, viz. r b = Kph 
vs. ti = qK e pc/(dv/dr) (see §2.2). Namely in both 
cases, the correction factor takes the form (1 — e~ r )/r; 
cf. eqns. © and lf3l)|) . 

Now, as noted in §2.2, in the theory of line-driven 
winds a crucial extensio n over early, single-line analyses 
(|Lucv fc Solomonlf l9701 was the introduction by CAK of 
a formalism for treating the cumulative effect of many 
lines. 

Likewise, recognizing the unrealistic nature of the 
above simple picture of a clumped medium in which all 
the blobs have identical properties, let us next consider 
a model in which the structured component consists of 
an ensemble of individual clumps with a range of optical 
depths Ti. Compared to a smooth, unclumped medium, 
the effective opacity of the clump component should then 
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be reduced by a factor 

h 



dT Tr 



1 - e 



(54) 

where the latter approximation assumes the summation 
can be approximated by an integral over a suitably de- 
fined clump distribution in optical depth, df / dr. 

In direct analogy with the CAK power-law formalism 
for line-strength, let us specifically make the ansatz that 
this distribution can be described by an exponentially 
truncated power law (see §§2.2 and 2.3, and eqn.|!||), 



1 



A _ 
dr r[a p ] \r, 



-t/t 



(55) 



where r now represents the optical depth of the strongest 
clump, and a p > is the power-index 5 . This assumed 
distribution obeys the normalization 



dr 



( 1 
dr 



1. 



(56) 



Applying eqn. 1551) in eqn. I|54|) . we find the opacity 
reduction of the clumped medium to be given by (cf. 
ean. lTU|l . for a p ^ 1, 



(1 



\l-a„ 



(1 - a p )r 
[(1 + 1/r,,) 1 ^ 



(1/To) 



1— ct v 



(1 - a p )To" 
or for the special case a p = 1, 



(57) 



(58) 



For fixed clump characteristics, the optical depth r D 
again (cf. eqn. l3*6"|) scales with the density p, 



Po 



(59) 



where the critical density at which this strongest blob 
has unit optical depth is (cf. eqn. [^7J| 



Po 



1 

nh ' 



(60) 



with h a the associated porosity length. This model of a 
medium with a power-law porosity is now characterized 
by two parameters, the power index a p , and the porosity 
length h Q of the strongest blob. 

In the limit that even the strongest clump is optically 
thin, t„ < 1, we recover from eqns. l|571) and that 
there is only a small reduction in opacity, 



h[T ] wl - a p r /2 ; t < 1. 



(61) 



In the opposite limit that the strongest clump is opti- 
cally thick, t 3> 1, the asymptotic scaling depends on 

5 We assume that ct p > because otherwise the contribution 
from wea k structure causes the overall number normalization inte- 
gral 1561 to diverge, unless one introduced a further free parameter 
to cut off the distribution at some minimum porosity length. The 
subscript "p" stands for porosity, and is added to distinguish this 
exponent from the usual CAK line-opacity power-index a. Finally, 
note that r[a p ] represents here the Gamma function, and not the 
Eddington parameter. 



whether a p is larger or smaller than one. For a p < 1, 
the square bracket term in f57)l approaches unity in this 
thick limit, and we obtain for the opacity reduction, 



h[To] 



1 



(1 - a p )r '• 
For a > 1, we find 



h[To] 



1 



(dtp - 1)t 



t > 1 & a p < 1 . (62) 



t > 1 & a p > 1. (63) 



As seen from eqn. (|58|) . the special case a p = 1 also scales 
as kb(r ) oc l/r in this optically thick limit, but with a 
proportionality factor that increases logarithmically with 
optical depth, i.e., ln(l + t g ). 

We thus see that the effective opacity of the structured 
component can become arbitrarily small at large densi- 
ties, p/ Po — t 3> 1. However, any "interclump" , smooth 
compoent should still contribute a constant opacity, in- 
dependent of density, so that the overall effective opacity 
is given by 

k[T ] = k min + k [T o ] , (64) 

which thus has a "floor" or minimum effective opacity, 
kmim set the mass fraction of any smooth, interclump 
component. 

6.2. Mass Loss Rate for the Power-Law Porosity Model 

Let us now apply this power-law porosity model to- 
ward computing the mass loss properties from a super- 
Eddington surface layer. In analogy with eqn. we 
first write the generic mass loss rate in a form that is 
scaled by the porosity-length h a and the associated crit- 
ical density p a = 1/k/i d , 



M = 



1 1 a*c 



(65) 



The sonic density is now found through application of 
eqn. 157|) or (|58|l to the sonic point condition, 



rk[p*/p ] 



(66) 



where now we must also require the mass fraction of in- 
terclump medium k m in < 1/T, so that eqn. I|66|l has 
a solution. In the following analyses, we assume this is 
always the case, regardless of how large T becomes. In 
effect, we thus henceforth effectively take 



k[p/p ] « h[p/p 



< 



P* 



(67) 



which applies from the sonic point outward, i.e., through 
any resulting wind outflow. 

In such a medium in which the unclumped Edding- 
ton parameter T is sufficiently above unity, this sonic 
point is located where the density is high enough that 
the strongest clump is quite optically thick, t S> 1. Un- 
der this circumstance, the structured component opacity 
reduction given by the full form 157|) can be reduced to 
the simple power-law forms l10'2[| and (|63|l , thereby allow- 
ing trivial solution of the critical condition l|66|) 



Po 



r 



1 - ou 



a p — 1 



l/a p 



; r > 1 k a p < 1 , (68) 
; r > 1 & a p > 1 . (69) 
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In the opposite limit for which the Eddington parameter 
is only marginally above unity, we find through applica- 
tion of eqn. (|ol")l in cqn. (H2jJ), 



p» _ 2(r- 1) 

Po Oi p 

For the special case a p — 
eqn. I|66[l is 

P 



K 1 . 



(70) 



1/2, a general solution to 



'— = 4r(r - 1) 

Po 



a p = 1/2 



(71) 



For the special case a p = 1, application of eqn. H58f) m 
l|66|) also yields a general, explicit solution, 



— = -i-rw_i[-exp(-i/r)/ri 

Po 



2riog(r) ; a p =l, 
(72) 

where W—i is now t he lower branch of the Produc t-Log 
(Lambert) function l|Jeffrev. Hare, fc Corliss! 1 1996(1 . and 
the latter approximation is roughly applicable for mod- 
erate r, e.g., remaining within ca. 25% of the correct 
solution for T < 10. 

Next, let us again introduce a porosity-length param- 
eter, rj = h /H, relating now to the strongest clump 
(cf. eqn. gU|). Applying this, plus eqns. (JHSJl-O, into 
cqn. i|65[) . we nn d the mass- loss rate for the power-law 
porosity model has the scalings 



M, 



1 



pow ' 



r 



1 - Ctr. 



l/a p 



L * 



Op - 1 

2(r-i) 



2io g (r) 
4(r - i) 



■n a*c 
r] a*c 
rj a*c 



T > 1 & a p < 1 (73) 
r > 1 & a p > 1 (74) 

r - 1 « i 



j 1 



a p = 1/2 



(75) 
(76) 
(77) 



It is of interest to compare these results with those 
given in eqn. Il l.'ill for the single-length porosity model 
of §5.2. For a p > 1, eqns. (|ZU and ((ZSJ) show that 
the scaling is similar to the single-length model in both 
the strongly (r » 1) and weakly (r — 1 <C 1) super- 
Eddington limits. As noted above, this stems from the 
fact that in such a p > 1 models the opacity reduction is 
dominated by structure with the largest porosity length. 

However, for a p < 1, comparison of eqn. (|73l) and 14311 
shows there can be a substantially larger mass loss in the 
strongly super-Eddington limit of the power-law model. 
In the single-length model, the opacity in dense layers is 
reduced by the inverse of the density, but in the a p < 
1 power-law models the reduction is weaker, scaling as 
l/p ap . For a large overall Eddington parameter T ^> 1, 
the sonic-point reduction of the effective Eddington to 
unity occurs at a deeper, denser layer, implying a larger 
mass loss rate. The net result is to increase the sensitivity 
of the mass loss rate to the Eddington parameter T. For 
example, for a p = 1/2, the overall scaling is in proportion 
to 4(r - 1) instead of 1 - 1/r 2 . 

Finally, let us consider this power-law-porosity mass- 
loss relative to the tiring limit discussed in §3. From the 



discussion in §5.2 for the single-porosity length model, 
eqn. (|47|l already gives the scaling of the tiring number 
associated with the basal mass loss rate, M* = L*/a*c. 
In the power-law models, the full tiring number is just 
increased by the additional dimcnsionless factors from 
eqns. Ij73(l - I|77|l . For example, for the analytic case a p = 
1/2, we find the tiring number scales as 



l pow ■ 



M pow 

i4 (r-i)^ 

Vo 



= 0.13 



r-i m^/r* 



1/2. (78) 



rj a 2 n M Q /R Q 
We therefore see that, under the canonical assumption 
that rj w 1 (i.e., h Q w H) and a2o ~ 1, a porosity model 
with power index a p = 1/2 and a moderately large Ed- 
dington parameter can readily approach the tiring limit, 
e.g., m pow ~ 0.52 for V = 5. 

Eqn. (j78(l suggests that while photon tiring should have 
marginal influence for only mildly super-Eddington mod- 
els, i.e., r < 2, it should limit the effective mass loss in 
strongly super-Eddington cases, i.e., T > 3. We examine 
this issue further in conjunction with determination of 
the associated wind velocity law (§6.4). 

6.3. Wind Velocity Laws for Power-Law Porosity 
Models without Tiring 

In addition to mass loss rate, the power-law porosity 
formalism can also be used to model the outward wind 
acceleration and resulting velocity law. Let us first exam- 
ine this under the assumption that the mass loss rate is 
small enough that photon tiring can be neglected. Ignor- 
ing as before gas pressure effects, the wind acceleration 
can now be described by the dimensionless equation of 
motion, 

w'(x) = Tk[T (x)] - 1. (79) 

As above, w' is the wind acceleration in units of the local 
gravitational acceleration, and x is a spatial coordinate 
defined by x = 1 — i?*/r, where r is the local radius and 
i?» is a characteristic wind sonic-point radius at which 
the RHS crosses zero, i.e., w'(0) = rfc[r o (0)] -1 = 0. 
Applying the opacity reduction eqn. I|57(l . the equation of 
motion l|79|) becomes an ordinary, first-order differential 
equation that can be integrated using standard numerical 
techniques. The free parameters in this integration are 
the power-law index a p , the Eddington parameter T, and 
the initial flow energy w(0) — u>*. 

Figure 2 shows results for the ratio of wind termi- 
nal speed to escape speed Voo/v esc (= y/ w[l]), plotted 
vs. Eddington parameter T for various a p and w*. The 
terminal speeds are typically of order the escape speed, 
with the larger values for larger T, larger a p , and lower 
w*. 

We also generally find that the height variation of ve- 
locity can be reasonably well fit by a modified "beta" 
velocity law of the form 



r{r) - ( r~; - n 2 )[l--y- 



r 



2;'» 



(80) 
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Fig. 2. — Ratio of terminal speed to escape speed, foo/fesc 
[= \J w(X) ], plotted as a function of Eddington parameter T, for 
selected values of ce p = 0.4, 0.5, and 0.6 (lower, middle, and upper 
curves of each triad), and for initial energy w„ = 0.01 (lower triad) 
and w* = 0.001 (upper triad). The results can be roughly fit by 
Voo/Vesc ~ [log(l/ TO *)(r - 1)/3]( 1 +«p)/ 4 . 




Fig. 3. — Effective velocity law index /3 e g (= 0.5 log[(to(l) — 
u>*)/(to(0.1) — to*)]) plotted as a function of Eddington parameter 
r, for selected values of a p = 0.4, 0.5, and 0.6 (lower, middle, and 
upper curves of each triad), and for initial energy to* = 0.01 (upper 
triad), to, = 0.003 (middle triad), and u>* = 0.001 (lower triad). 



with however the velocity power index (3 dependent on 
the parameters T, a p , and td*. To characterize this de- 
pendence, we define an effective index 



(3 cS = 0.5 log 



w(l) 



iu(0.1) 



(81) 



which we find gives a reasonable fit to the numerically 
integrated velocity variation. Fig. 3 plots f3 c s vs. T for 
various a p and u>*. The values are typically within 25% 
of unity, with larger /3 e fr occurring for larger T, larger u>*, 
and smaller a p . 

By comparison, for line-driven winds typical terminal 
speeds are 2-3 times the escape speed, with a velocity- 
law power index /3 0.8 — 1. Figs. 2 and 3 show similar 
results for the power-law porosity model of continuum 
driving, with however a somewhat wider range that de- 
pends on the parameters T, a p , and ttj*. 

6.4. Photon Tiring and Wind Stagnation for Strongly 
Super-Eddington Porosity Models 
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Fig. 4. — a.) Velocity over escape velocity plotted vs. scaled 
inverse radius x, for Eddington parameters T=4 and 5, for models 
ignoring photon tiring (m* = 0) and with basal tiring parameter 
m» = t)g SC /2a»c = 0.032. b.) Terminal speed over escape speed, 
plotted vs. Eddington parameter T, for models with basal tiring 
parameter m» = 0.032, with various power-indexes a p and initial 
energies to*. 



As noted in §6.2, the mass loss rates expected for 
power-law porosity models with moderately large Ed- 
dington parameters (T > 3 — 4) approach the photon tir- 
ing limit. Let us now derive velocity solutions for poros- 
ity models that take into account photon tiring. From 
eqns. i|27|l and (|57J) the equation of motion now takes the 
form 

w'(x) = -1 + k[r (x)] r [1 - m(w + x)]. (82) 

For a given tiring number m, this can again be inte- 
grated numerically from the initial value u>», given also 
the model parameters T and a p , 

Using eqns. 147|) and 165f) . we can write the tiring num- 
ber as 

m = (83) 

where m* = v^ sc /2a.^c — a«/2w* c, a nd p*/p is a set by 
a and T through solution of eqn. I|66[) . For the typical pa- 
rameters a p — 1/2 and ui* = 0.001, figure 4a then plots 
the velocity over escape velocity, v/v eS c, vs. scaled in- 
verse radius x, comparing results without tiring (to* = 0) 
and with the fiducial value of the basal tiring number 
777* = 0.032 (ean. I47fl . Note that without tiring the ve- 
locity is higher for the higher Eddington parameter, but 
with tiring the velocity is reduced, with the stronger re- 
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Fig. 5. — Terminal speed over escape speed, plotted vs. photon 
tiring number m = M / Mu r , for models with various specific power 
indices over the range a v = 0.4 — 0.9 (in increments of 0.1). The 
dashed curves show contours of the associated Eddington parame- 
ter, for the specific values T = 2, 3, 4, 5, and 10. The convergence 
of models along the upper right bound represents the total tiring 
limit. 



duction now making the velocity lower for T = 5 than 
for T = 4. 

For this same basal tiring number m* = 0.032, Fig. 4b 
shows the terminal speed over escape speed, Voo/ W esc, 
plotted vs. the Eddington parameter L, for selected pa- 
rameters a p and w*. Again the tiring reduces the speed, 
but note that for the lower power index, a p = 0.4 with 
moderately large T > 3.5, the wind can no longer reach 
large radii (x = 1) with a finite speed. 

Fig. 5 plots the ratio of terminal speed to escape speed 
vs. photon tiring number m = M/Mu r , for models fixed 
with the standard parameter set r\ = 1, m» = 0.032, and 
= 0.001, but with various specific power indices over 
the range a p = 0.4 — 0.9 (in increments of 0.1). The 
dashed curves show contours of the associated Edding- 
ton parameter, for the specific values L = 2, 3, 4, 5, and 
10. Note that for lower a p and moderately large L the 
tiring becomes substantial. Indeed, the convergence of 
models along the upper right bound represents the total 
tiring limit, for which the combination of kinetic and po- 
tential energy of the wind has exhausted the entire base 
luminosity L». Hence, models along this limit have little 
or no remaining luminosity to be observable as radiation. 

This illustrates a remarkable property of these power- 
law porosity models. For small power exponents and 
moderately large Eddington parameters, they can readily 
exhaust nearly the entire available luminosity in driving 
the mass loss. 

7. DISCUSSION 

7.1. Comparison with Previous Porosity Analyses 

Let us now consider how the above scalings com- 
pare with results from previous analyses of porosity- 
moderate d mass loss. S pecifically, as noted above (cf. 
can. I45fl . Shaviv (2001bJ) has cast the mass loss scaling 
in terms of a "wind function" W(r) times the deviation 
of the luminosity from its Eddington limit value, 

M = W[T] L * - LeM = W[T] fc^) — • (84) 



As noted in §5.2, for the single-scale porosity model, the 
wind function W[T] « (T + 1)/Trj, which therefore varies 
from VV ~ 2/rj for the marginally super-Eddington case 
r-l«ltoW«l/i| for t he strong super -Eddington 
limit r ^> 1. For comparison. IShavrvl l)2001bD considered 
two limiting cases, of either blobs in the optically thin 
limit or thick but vertically elongated, and a third case 
of a two phase Markovian mixture. For all cases, he 
estimated that the wind function would be of order unity, 
which is indeed our estimate here from the single-scale 
porosity model. 

However, for porosity models with power-indices a < 
1, we find here that there can be a further increase above 
this basal mass loss of the single-scale model. For exam- 
ple, for the canonical case a p = 1/2, comparison of l|84|l 
and (|77J) implies the wind function scaling to be 

4r 

W[L] = — ; op = 1/2. (85) 
r\o 

A key new feature of the power-law porosity model is 
consequently that this wind function increases with the 
Eddington parameter, and so now can become substan- 
tially larger than just the order unity value derived for 
the sin gle-scale model, and also inferred from the anal- 
yses bv lShavivl l|2001bD . Since even the basal mass loss 
rate l|39l) for the single-model can be a few percent of 
the tiring limit (|20|l . the additional increase for moder- 
ately large Eddington parameter T > 4 — 5 can bring the 
overall mass loss close to this limit. 

In addition, of course, the mass loss scaling depends 
on the scaling of the porosity- length h a , and this might 
also depend on T. In particular, we have assumed above 
that the porosity length scales with the pure-gas-pressure 
scale height H] but in a medium in which there is a 
residual smooth component that puts a floor k m i n in 
the effective opacity (see §§5.3-6.2), the effective grav- 
itational scale height at large depth is instead given by 
H/(l — kminX)- Scaling the porosity length scales with 
this larger stratification height could then reduce the 
mass loss rate by a factor (1 — k m i n T), which in the 
possible case that k m i n T < 1, could represent a sub- 
st antial reduction . This general effect was also discussed 
bv lShavivl <l2001bfl . 

7.2. Comparison to Inferred Mass Loss of Giant 
Eruptions 

Let us next consider the conditions needed for this 
power-law-porosity model to reach the observationally 
inferred mass loss rates of giant eruptions of Luminous 
Blue Variable (LBV) stars like 77 Carinae. A nalyses of the 
resulting Homunculus nebula l)Smitbll2002ft indicate that 
during the roughly 20-year outburst from ca. 1840-60, 
r\ Carinae lost 2 — 10M©, representing an average mass 
loss rate of 0.1 — 0.5M Q /yr. Moreover, the expansion 
of the Homunculus nebula is inferred to be in the range 
Vexp = 500 — 700 km/s. If we characterize this expansion 
velocity as being some order-unity factor of the escape 
speed from the surface of origin, ^Woo = v exp /v e sc, then 
the total luminosity associated with this wind mass loss 
is 

Lwind = Mv 2 esc (1 + Woo )/2 

= 30 x 10 6 L o Mi ™; IR * (l + woc), (86) 
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Fig. 6. — Same models as fig. 5, but now for corresponding di- 
mensional values, showing terminal speed vs. mass loss rate (scaled 
by Mioo, the stellar mass in units of 100 Mq). The dashed curves 
again show contours of the Eddington parameter, now for integer 
values in the range V = 2 - 20. Both plot sets assume the same 
dimensionless parameter combination, w* = 0.001 and ra* 0.032, 
which together imply dimensional parameters of a* = 20 km/s 
for the base sound speed and v esc = 620 km/s (= u esc ,0) for the 
surface escape speed. The semi-transparent box denotes the obser- 
vationally inferred parameter range for the 1840-60 giant eruption 
of r\ Carinae that gave rise to the Homunculus nebula. 



where M% = M/(lM Q /yr). For a solar value of mass to 
radius, v esc ~ 620 km/s pa v exp ; thus, assuming pa 1 
in this case, we see from eqn. (|86[) that the observation- 
ally inferred range of mass loss for r\ Carinae requires a 
total wind luminosity in the range L W i n d pa 6— 30 x lQ e Lq 
ljSmithM2002D . 

Moreover historical observations suggest a radia- 
tive luminosity of roughly L b a ~ 20 x 1O 6 L 
ijDavidson &: Humphreys 111997(1 . implying a total lumi- 
nosity of up to 50 x 10 6 L Q ! 

To facilitate comparisons with these observational val- 
ues, let us now convert the above dimensionless model 
results into physical units. The above standard di- 
mensionless parameters w* = (u esc /a*) 2 = 0.001 and 
m* = Vg SC /2a*c = 0.032 imply specific dimensional val- 
ues for the surface escape speed 

v esc = Icm^^fwl pa 620km/s pa v esc , Q , (87) 

and the surface sound speed 

a* = 2cm»u;* pa 20km/s. (88) 

The dimensional value of the terminal speed is thus given 
by Voo = 620 km/s \J From eqn. I|38|) . the dimen- 
sional mass loss is obtained from 

M = 3 .9xl0- 3 ^^^, (89) 
yr p rj o ci20 

where M 100 = M*/100M Q and a 20 = a*/(20 km/s). 

Fig. 6 plots the terminal speed vs. the mass loss 
rate M for the same models shown in the dimensionless 
plots of fig. 5, for various power indices a p . The dashed 
curves again show contours of the associated Eddington 
parameter T. The observationally inferred ranges of mass 
loss and terminal speed are denoted by the shaded box. 
Assuming a canonical stellar mass Mioo = 1, note that 
even the upper range of mass loss, M\ pa 0.5, can be 
obtained with a v ~ 0.6 and T pa 8.5. 



However, note also from eqn. (|86(l that the associated 
wind luminosity for this case is L W i n d pa 30 x 10 6 L Q . 
From the conversion T pa Q.27Lq/Miqq, this implies an 
associated "wind Eddington parameter", T W i nc i sa 8.1, 
which therefore means that a fraction 8.1/8.5, or about 
95%, of the source luminosity is consumed in the wind! 
This in turn implies that the observable radiative lumi- 
nosity in such a model would be a mere 

L r ad, oc ~ (T ^wind)-^Edd 

pa (8.5 - 8.1) 3.7 x 10 6 L© 

«1.5xlO 6 i , (90) 

which is much less than the historically ob- 
served value of L„b R ~ 20 x 10 6 L© noted above 
(pavidson fc Humphreys lll997Tl . 

Our further analysis indicates that fitting this observed 
luminosity together with the inferred mass loss and ter- 
minal speed requires a reduction of the assumed surface 
escape speed. As shown in fig. 2, without photon tiring 
the ratio of terminal speed to escape speed increases with 
increasing Eddington parameter T. This is a direct re- 
sult of the fact that radiative driving scales with T, and 
so as seen from the equation of motion without tiring 
(eqn. I79fl the scaled acceleration likewise increases with 
r. But the mass loss also increases with T, and in models 
with photon tiring, the increased tiring number reduces 
the driving (see ean.l82|l. Thus, as shown in figs. 4 and 5, 
the terminal speed peaks and then decreases with higher 
T. This makes it possible to obtain both a large mass 
rate and a modest terminal speed that is equal to the as- 
sumed escape speed. But because of the dominant role of 
tiring, a by-product of such models is that there is little 
remaining radiative energy to be observed as emergent 
luminosity. 

However, if we assume a lower escape speed, then 
models with less tiring can produce the observed ter- 
minal speed. For example, taking the same stellar mass 
-Mioo = 1 but an increased stellar radius i?ioo = 3 im- 
plies a reduced effective escape speed v esc = 360 km/s, 
so that the same target expansion speed now requires a 
speed ratio foo/^esc = V% ~ 1.73, and therefore an en- 
ergy ratio Woo = 3. For the targeted mass loss Mi = 0.5, 
we then see from eqn. that the associated wind lu- 

minosity is now L w i n d ~ 20 x 10 6 Lq. Together with the 
observed luminosity, this requires of total base luminos- 
ity of L* pa 40 x 10 6 -Lq, corresponding to an Eddington 
parameter T pa 10.8. 

By taking this with Mioo = 1 and i?ioo = 3, 
and keeping the remaining three parameters the same 
(rj = 1, a = 0.6, a2o = 1), we then find that the ob- 
servable radiative luminosity is now higher than before 
(L bs ~ 13.4 x 10 6 L©), but still below the target value, 
essentially because of too much tiring from the still-too- 
high mass loss rate (Mi pa 0.76). To reduce the mass 
loss, we can increase the porosity parameter r\ Q and/or 
the sound speed 020- Choosing the former, we find that 
rj = 1.56 gives results within 2% of the targeted values 
for all three observational parameters. 

Of course, given that the model has a total of six ad- 
justable parameters, the significance of being able to 
fit three observed values should not be exaggerated. 
Nonetheless, it is encouraging that this power-law poros- 
ity model can reproduce such rather extreme mass loss 
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conditions with quite plausible choices of these param- 
eters. It is particularly interesting that, because of the 
scalings of mass loss and flow speed with the Edding- 
ton parameter, and their interplay with photon tiring, 
reproducing both the inferred mass loss and radiative lu- 
minosity requires that the source star's surface escape 
speed is somewhat smaller than the inferred expansion 
speed. 

7.3. Gravity Darkening and the Shaping of LBV 
Nebulae 

For simplicity, the above formalism has assumed a 
spherically symmetric mass loss, but in the case of 
77 Carinae, the resulting Homunculus nebula exhibits 
a distinctly bipolar, prolate form. Detailed analysis 
( Smith] 12002(1 indicates the nebula has both a higher 
speed and greater density toward the bipolar symme- 
try axis, suggesting that the 1840-60 giant outburst 
was itself intrinsically bipolar. Moreover, there is evi- 
dence, h^HL^StiLSiSliSi^: 86 ^ intcrfcrometric observa- 
tions i|va n Boekel et al]l2003() and from HST slit-spectra 
( Smit h et al]l2003lK the present-day wind of 77 Carinae is 
also bipolar about the same symmetry axis as the Ho- 
munculus. 

One promisin g sce nario, first proposed by 
Owocki fc Gavlevl lll997l) and developed furt her by 
Maeder fe Mevnetl (12000). iMaeder fc Desiacauesl (|2000|) . 
and IDwarkadas fc Owockil |J2002|) . is that such bipo- 
lar mass loss is the natural consequence of radiative 
driving from a nearly critically rotating s tar with a 
subst antial equatorial gravity darkening ijvon Zeipell 
1924). However, since this basic idea arose from 
models of relatively low-density , line-driven winds 
l|Owocki. Cranmer. fc Gavievlll996|) . the applicability to 
much denser, continuum-driven models for r\ Carinae or 
other LBVs has been uncertain. Eventually, the issue 
should be examined through 2-D dynamical models that 
account for the latitudinal mass and radiation transport 
within the optically thick, bipolar expansion. 

But for now, it is worth noting that the underlying 
scalings for the speed and mass flux have some key sim- 
ilarities in both the basic line-driven-wind theory and 
the present models for porosity-moderated continuum- 
driving of a dense, optically thick flow. In particular, 
comparison of cqns. l|16fl and (JTSJ-JTTJ) shows that the 
mass loss rate in both cases scales with the stellar lumi- 
nosity times a correction factor that is a function of the 
Eddington parameter, i.e., M cx L*/[r]. Accordingly, 
in a rotating star, the local surface mass flux m*{9) at 
any colatitude 9 should vary in proportion to the local 
surface radiation flux F*(6), times a function of the lo- 
cal effective Eddington parameter T = F^/g c s, where 
<7eff is the effective, c entrifugally reduc ed surface gravity. 
But for the standard lvon Zeipell l|1924f) gravity darkening 
scaling that F*(6) cx g e s(0), we see that the Eddington 
parameter is latitudinally constant, implying then that 
the surface mass flux should scale simply as 

m*(0) ocF* (0) ocg eS (8). (91) 

Since both the radiative flux and effective gravity are 
maximum at the rotational pole, eqn. (j91(l shows that 
the mass flux should be strongest near the poles. 

Similar arguments can be made for the latitudinal vari- 
ation of the flow expansion speed. Again, the detailed re- 



sults may depend on latitudinal components of the mass 
flow, radiation flux, and radiative force, and should be 
eventually be analyzed through 2-D models. But within 
the context of simple 1-D scaling relations for both line- 
driven and porosity models, the outflow speed should 
follow the approximate scaling, 

Uoo (0) cx u eac (0) cx y/g cS (0). (92) 

Since the effective gravity is highest toward the poles, 
we can expect the nebula expansion to be faster near the 
symmetry axis. Observations do indeed show that v{9) 
for the polar lobes of the Homunculus nebula is roughly 
proportional to the simple l atitudinal va riation of escape 
speed from a rotating star (|Smithll2002|) . 

Overall, the expected faster polar flow speed can ex- 
plain the generally prolate form of the expanding neb- 
ula, while the higher polar mass flux can explain the 
observationally inferred mass concentration near the po- 
lar symmetry axis. Thus, an attractive feature of the 
porosity-moderated continuum-driven formalism is that 
it preserves these key 1-D flow scalings from line-driven 
models, while allowing extension to much larger mass 
loss rates. However, as noted, more complete 2-D mod- 
els should be developed to examine how these general 
scalings might be affected by latitudinal mass and radi- 
ation transport. 

Finally, while several LBV nebulae show such a bipolar 
or ellipsoidal form, none except maybe HR Carinae is as 
severely pinched at the equator as 77 Carinae. P Cygni 

- the only other Galactic star observed to have a giant 
eruption approaching the extremity of th at in 77 Carinae 

- has a nearly spherical ring nebula ijNota fc Clampir] 
119971) . In this context, it should be emphasized that the 
basic porosity model developed here can accomodate a 
range of geometric forms, depending on the degree of 
rotation of the source star. 

8. CONCLUDING SUMMARY 

Let us conclude with an itemized summary of the goals, 
methods, and results of the analyses in this paper. 

1. The overall goal is to develop a radiative-driving 
formalism for explaining the extremely large mass 
loss rates and moderately high outflow speeds in- 
ferred in giant outbursts of Luminous Blue Vari- 
ables, most particularly the 1840-60 outburst of 
77 Carinae that resulted in the Homunculus neb- 
ula. The general approach combines and builds 
upon two previous models for radiative driving, 
namely the well-established CAK line-driving gen- 
erally applied for more quiescent phases of hot-star 
wind mass loss, and recent notions of porosity- 
moderated continuum drivin g in stars that exceed a 
generalized Eddington limit (|shavivlll"998l I2001b|) . 

2. As a basis for this synthesis, we first review the 
CAK line-driven formalism, showing thereby that 
for reasonable values of the line-opacity normaliza- 
tion, the expected wind momentum (i.e., the prod- 
uct of mass loss rate and flow speed) in such models 
falls well below what is inferred in LBV giant out- 
bursts. The formal divergence in CAK mass loss 
rate as a star approaches the Eddington limit is 
accompanied by a vanishing terminal speed, and 
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moreover in practice is limited by the "photon tir- 
ing" of the finite luminosity available to lift mate- 
rial from the gravitational potential at the stellar 
surface. 

3. In the deep stellar interior, approaching and ex- 
ceeding the Eddington limit leads to convective en- 
ergy transport that lowers the radiative flux and 
allows a normal hydrostatic stratification with out- 
wardly declining pressure and density. In near- 
surface layers where the lower density makes con- 
vective transport inefficient, a super-Eddington 
condition again implies that the outward radiative 
force exceeds gravity. Hydrostatic stratification 
then requires an outward inversion of pressure and 
density, which however can only be maintained over 
a limited range, given the inherent outer boundary 
condition of vanishing pressure and density. 

4. Moreover, any outflow initiated where convection 
first becomes inefficient would imply a huge mass 
loss rate, of order L»/a*, where a* is the sound 
speed; this exceeds by a large factor the "tiring 
limit" mass loss, Mtir — 2£*/t>e SC , for which the 
energy expended to lift material from the surface 
gravitational potential equals the assumed stellar 
luminosity. The flow stagnation associated with 
photon tiring, together with other instabilities, 
seems likely to impart a complex, time-dependent, 
3-D structure to the atmosphere of any star that 
approaches or exceeds the Eddington limit. 

5. As first noted bv lShavivl l)1998|) . the "porosity" of 
such a structured, super-Eddington medium can 
lower the effective driving opacity in the deeper, 
denser layers, allowing a base hydrostatic balance 
that transitions to a supersonic outflow as the 
structures become optically thin, and hence are 
again subject to the full radiative force. 

6. We introduce a simple formalism that character- 
izes this effect in terms of a "porosity length" h = 
L 3 /l 2 , where I and L represent the size and separa- 
tion of individual clumps or blobs. We show that a 
simple ansatz (somewhat analogous to the mixing 
length parameterization of convective transport) - 
that this porosity length is an order-unity factor n 
times the gravitational scale height H - leads to a 
mass loss rate M por w (1 — l/T^L^/na^c, a scal- 
in g that is very s imilar to that derived previously 
bv IShavivl (|200lUl for both vertical elongation and 
Markovian mixture models of atmospheric struc- 
ture. Though large, such mass loss is still typically 
only a few percent of the tiring limit. 

7. In analogy with the CAK formalism of driving by a 
power-law ensemble of lines, we then generalize to a 
"power-law-porosity" treatment, characterized by 
a power index a p and a porosity-parameter r\ for 
the most optically thick blob. For a p > 1, the de- 
rived mass loss rates are similar to the single-scale 
model, but for a p < 1 they lead to enhancement by 
factors that scale with the Eddington parameter as 
p-i+i/ap^ -p Qr i ar g e Y 7 small dp, and/or small rj , 
this can now lead to overall mass loss rates that 
approach the photon tiring limit. 



8. Without tiring, the wind velocity follows a canon- 
ical "beta- velocity-law" form, with index j3 « 1, 
and a terminal speed proportional to the escape 
speed times a factor that scales roughly with (T — 
l)( Q p +1 )/ 4 . With tiring, the velocity law can be- 
come nonmonotonic, with a terminal speed that 
decreases with increasing mass loss rate. Models 
with a terminal speed less than or equal to the sur- 
face escape speed are strongly tired, and so have a 
greatly reduced observable radiative luminosity. 



9. Given stellar parameters L„, M*, and plus a 
base sound speed a*, the mass loss properties of 
the power-law-porosity model depend on the power 
index a p and porosity-length parameter n . For 
quite reasonable values, the model can reproduce 
the observationally inferred mass loss and speed 
for the giant eruption of n Carinae, but matching 
also the historically estimated radiative luminosity 
during this epoch requires a modest surface escape 
speed, with ratio of stellar mass to radius perhaps 
a third of the solar value. 



10. Though the magnitude of mass loss greatly ex- 
ceeds what's feasible with line-driving, the poros- 
ity model retains the key scalings with gravity and 
radiative flux that would give a rapidly rotating, 
gravity-darkened source star an enhanced polar 
mass loss and flow speed, as is inferred from the 
bipolar form of the Homunculus nebula in n Cari- 
nae. 

Overall, the extended porosity formalism developed 
here provides a promising basis for self-consistent dynam- 
ical modeling of even the most extreme mass loss out- 
bursts of Luminous Blue Variables, namely those that, 
like the giant eruption of n Carinae, approach the pho- 
ton tiring limit. But of course much further work will be 
needed to substantiate and quantify the basic "porosity- 
length" phenomenology proposed here. Much in the 
same way that modern-day multi-dimensional simula- 
tions of stellar convection have been used to test clas- 
sical mixing length phenomenology, future work should 
focus on developing 2-D or 3-D radiation hydrodynami- 
cal simulations of the nonlinear evolution of atmospheric 
structure in stars near and above the Eddington limit. A 
particular challenge will be to develop fast techniques for 
treating the multi-dimensional, nonlocal radiation trans- 
port in such a medium. Of course, until such more fun- 
damental calculations are carried out, it is difficult to ap- 
praise the overall applicability of our porosity-length ap- 
proach. But in any case the basic phenomenology helps 
provide a motivation and conceptual framework for car- 
rying out and interpreting such more fundamental and 
challenging radiation hydrodynamical simulations. 
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